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Abstract. Of concern in the paper is generalized a theoretical study concerning the peristaltic 
flow of blood in the micro-circulatory system. The vessel is considered to be of non- uniform cross- 
section and blood to be a non-Newtonian fluid. The progressive wave front of the peristaltic flow 
is supposed sinusoidal/straight section dominated (SSD) (expansion/contraction type); Reynolds 
number is considered to be small with reference to the flow of physiological fluids. The non- 
Newtonian behaviour of blood is illustrated by considering the Herschel-Bulkley fluid model. 
The objective of the study has been to examine the effect of the effects of amplitude ratio, 
mean pressure gradient, yield stress and the power law index on the velocity distribution, wall 
shear stress, streamline pattern and trapping. Considerable quantitative differences between 
the results obtained for transport in two dimensional channel and an axisymmctric circular tube 
are noticed. The study shows that peristaltic pumping, flow velocity and wall shear stress are 
significantly affected by the non-uniform geometry of blood vessels. Moreover, the magnitude 
of the amplitude ratio and the value of the fluid index are important parameters that affect 
the flow behaviour. The novel feature due to SSD wave on said flow behaviour has also been 
discussed. 
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1. Introduction 

Peristaltic transport through different vessels of physiological systems by means of peristaltic 
waves is well known to physiologists as a natural mechanism of pumping materials in the 
case of most physiological fluids. Apart from physiological studies, many of the essential 
fluid mechanical characteristics of peristalsis have also been elucidated in different engineering 
analyses carried out by different researchers. Such types of characteristics usually occur when 
the flow is induced by a progressive wave of area contraction/expansion along the length of 
the boundary of a fluid-filled distensible tube. Benefits that can be derived from regarding the 
studies on peristaltic movement are available in our earlier communications (Misra and Pandey 
1994, 1995, 1999, 2001a, 2001b, 2002, Misra et al 2008, maiti and misra 2011) and also in other 
references e.g., Guyton and Hall (2006), Jaffrin and Shapiro (1971). 



Nomenclature 




b 


Wave amplitude 


do 


Radius of the tube at the inlet 


H 


Displacement of the wall in the radial direction 


k 


Reciprocal of n 


h 


A parametric constant 


n 


Flow index number 


P 


Fluid pressure 


Q 


Flux at axial location 


h 


Time 


R, 6, Z 


Cylindrical co-ordinates 


u,v,w 


Velocity components in Z, R, 9 directions respectively 


5 


Wave number 


Ap 


Pressure difference between ends of the tube 


A 


Wave length of the travelling wave motion at the wall 


n 


Blood viscosity 


V 


Kinematic viscosity of blood 





Amplitude ratio 


P 


Density of blood 


TO 


Yield stress of blood 


T h 


Wall shear stress 



The phenomenon of peristalsis has quite important applications in the design and construction 
of many useful devices of biomedical engineering and technology, artificial blood devices such 
as heart-lung machine, blood pump machine, dialysis machine etc. A theoretical foundation of 
peristaltic transport primarily with inertia-free Newtonian flows driven by sinusoidal transverse 
waves of small amplitude was suggested by Fung and Yih (1968). Shapiro et al (1969) 
first presented the closed form solutions for an infinite train of peristaltic waves for small 
Reynolds number flow, where the wave length is long but of arbitrary wave amplitude. The 
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said investigations (Fung and Yih 1968, Shapiro et al 1969) witnessed a variety of important 
applications that have been proved to be quite useful in elucidating the functioning of various 
physiological systems, such as flows in the ureter, the gastro- intestinal tract, the small blood 
vessels of the micro-circulatory system and other glandular ducts. Conditions for the the 
existence of physiologically significant phenomena of trapping and reflux in peristaltic transport 
were also determined by Shapiro et al (1969). Some of the references to earlier literatures on 
peristaltic flow of physiological fluids are available in (Jaffrin and Shapiro 1971, Srivastava and 
Srivatava 1984), while brief reviews of relatively recent literatures can be found in (Tsiklauri and 
Beresnev 2001, Mishra and Rao 2005, Yaniv et al 2009, Jimenez- Lozano et al 2009). Several 
other studies concerned with the analysis of different problems of peristaltic transport of various 
physiological fluids were reported by different investigators (Usha and Rao 1995, Takabatake 
and Ayukawa 1982, Jimenez-Lozano and Sen 2010, Bohme and Friedrich 1983, Srivastava and 
Srivastava 1985, Provost and Schwarz 1994, Chakraborty 2006). Past experimental observations 
indicate that the non- Newtonian behaviour of whole blood mainly owes its origin to the presence 
of erythrocytes. As early as in the sixth decade of the last century, different groups of scientists 
(Rand et al 1964, Bugliarello et al 1965, Chien et al 1965) made an observation that the non- 
Newtonian character of blood is prominent as soon as the hematocrit rises above 20%. However, 
it plays a dominant role, when hematocrit level lies between 40% and 70%. 

It is known that the nature of blood flow in small vessels (radius <0.1 mm) and at low shear 
rate (< 20sec~ 1 ) can be represented by a power law fluid (Charm and Kurland 1965, 1974). As 
reported by Merrill et al (1965), Casson model seems to be satisfactory for blood flowing in 
tubes of 130-1000/i. Scott Blair and Spanner (1974) estimated that blood obeys Casson model 
for moderate shear rate flows. However, they did not report any difference between Casson and 
Herschel-Bulkley plots over the range where Casson plot is valid for blood. They observed that 
for cow's blood, Herschel-Bulkley model is more appropriate than Casson model has been found 
by them. The results presented in the sequel for shear thinning and shear thickening fluids are 
quite relevant for the study of blood rheology. It is important to mention here that normal 
blood usually behaves like a shear thinning fluid for which with the increase in shear rate, the 
viscosity decreases. As mentioned in (Fung 1981, White 1974, Xue 2005) in the case of hardened 
red blood cell suspensions, the fluid behaviour is that of a shear thickening fluid for which with 
the fluid viscosity is enhanced due to an elevation of the shear rate. 

It is well-known that physiological vessels are usually non-uniform (Wiedman 1963, 
Wiederhielm 1967, Lee and Fung 1971). Some initial attempts to perform theoretical studies 
pertaining to peristaltic transport of physiological fluids in vessels of non-uniform cross section 
were made in several researchers (Srivastava and Srivastava 1984, Wiederhielm 1967, Lee and 
Fung 1971, Gupta and Seshadri 1976, Srivastava and Srivastava 1983). In all these studies, the 
physiological fluids were considered either as Newtonian fluids or non-Newtonian fluids described 
as Casson/power-law fluids. However, the different flow characteristics have not been adequately 
explained in these studies. Herschel-Bulkley models are of course, more suitable to represent 
physiological fluids, because the fluids represented by this model describe very well material 
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flows with a non-linear stress relationship depicting the behaviour of shear-thinning/shear- 
thickening fluids that are of much importance in the field of biomedical engineering (Malek 
et al 2002). Among the various types of non- Newtonian models used to represent physiological 
fluids, Herschel-Bulkley fluid model is more general. Use of this model has the advantage that 
results for fluids represented by Bingham plastic model, power law model, Newtonian fluid 
model can be derived from those of the Herschel-Bulkley fluid model. Also Herschel-Bulkley 
fluid model is found to yield more accurate results than many other non-Newtonian models. 

In the consideration of all the above, a study concerning the peristaltic flow of blood in a 
non-uniform vessel of micro-circulatory system has been undertaken here, by treating blood as 
a non-Newtonian fluid. It is worthwhile to mention that although flow through axisymmetric 
tubes is qualitatively almost similar to the case of flow in channels, magnitude of different 
physiological quantities associated with flow and pressure differ appreciably differ by one another. 
The mathematical analysis presented here is particularly suitable for investigating the peristaltic 
motion of blood through vessels in micro- vessels, e.g. arterioles and venules. Since for blood flow 
in the micro-circulatory system, the Reynolds number is low and since the ratio between radius of 
the tube and the wave length is small, the theoretical analysis of the problem under consideration 
has been performed by using the lubrication theory (Shapiro et al 1969). Basing upon the study, 
extensive numerical calculations have been made with the help of derived analytical expressions. 
Keeping a specific situation of micro-circulation in view, with the purpose of examining the 
pumping performance, the computational results are presented for the velocity distribution of 
blood, wall shear stress distribution as well as the streamline pattern and trapping. Influence of 
SSD wave front on various flow behaviour of peristaltic movement has also been discussed. The 
plots give a clear view of the qualitative variation of various fluid mechanical parameters. 

The analysis bears a strong promise of some important applications by which we can have 
a better insight of the micro-circulatory system. The analysis has an important bearing on 
the clinical procedure of extra-corporeal circulation of blood through the use of the heart-lung 
machine that may result in damage of erythrocytes owing to high variation of the wall shear 
stress. In addition, the results should useful application in the development of roller pumps and 
arthro pumps by which fluids are transported in living organs. 

2. Formulation 

Peristaltic motion of blood in arterioles and venules will be studied here, by considering blood as 
an incompressible viscous non-Newtonian fluid. The non-Newtonian viscous behaviour has be 
represented by Herschel Bulkley model. Non-uniform cross-sectional nature of the arteriole and 
venule have been emphasized in the study in the study. However, arteriole and venule have been 
assumed to be axisymmetric and the flow that take place through it has also been considered 
axi-symmetric. 

We take (R,0,Z) as the cylindrical coordinates of the location of a fluid particle, R being the 
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Fig. 1 : A physical sketch of the problem 



radius of the tube. Z is measured in the direction of wave propagation. Let R = H (cf. Fig. 
1) represent the wall of arteriole or venule that is considered to be induced by a progressive 
sinusoidal, SSD expansion or contraction wave train respectively propagating with a constant 
speed c travelling down the wall, such that 



H 



H 



H 



2tt 

a(Z) + b sin( — [Z — cti)), 
A 

a(Z) + b sin(^-(Z-c*i) 
a(Z) 

a(Z) - b sin(^(Z - cti) 
a(Z) 



if 0<Z<\ C1 

if A Cl < Z < A, 

if 0<Z<X C1 

if A C1 < Z < A; 



(1) 
(2) 

(3) 



with a(Z)=a + k{L where a(Z) represents the radius of the arteriole or venule at any axial 
distance Z from inlet, do being the radius at the inlet, k\(<l) a constant whose magnitude 
depends on the length of the arteriole or venule as well as on the inlet and exit dimensions, b 
the wave amplitude, t\ the time and A the wave length. The SSD expansion/contraction wave 
defined by equations (2) and (3) means that these waves are confined to a portion of length A Cl . 

3. Analysis 



Withe the consideration made above, the motion of blood in the arteriole or venule can be 
assumed to be governed by the partial differential field equations 

V • V = (4) 



„ £ dV 

Va + pf = p—, 
at i 



(5) 



where V is the velocity, f the body force per unit mass, p the density and the material time 



derivative. Further, a being the Cauchy stress defined by 

a = —PI + T, 
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T = 2fiD + S, 
S = 2rjD, 

D being the symmetric part of the velocity gradient, 

D = \[L + L T lh= VV, 
—PI denotes the indeterminate part of the stress due to the constraint of incompressibility, 
where \i and rj represent viscosity parameters. The Herschel-Bulkley model is the representative 
of the combined effect of Bingham plastic and power-law behavior of the fluid. When strain-rate 
is low such that 7 < the material acts like a viscous fluid with constant viscosity /i . But 
when the strain rate increases and the yield stress threshold, f , is reached, the fluid behavior 
is described by a power law of the form 

'/ . * 

7 

where m and n denote respectively the consistency factors and the power law index, n < 1 
corresponds to a shear thinning fluid, while for a shear thickening fluid n > 1. In the case of 
a uniform tube, if the tube length is an integral multiple of the wavelength and the pressure 
difference between the ends of the tube is a constant, the flow is steady in the wave frame. Since 
in the present study, the geometry of wall surface is non-uniform, the flow is inherently unsteady 
in the laboratory frame as well as in the wave frame of reference. The equations (Herschel and 
Bulkley 1926) that govern the fluid motion in a micro-vessel : 

fdU dU dU\ OP 1 d(Rr RZ ) dr zz 

(dV TT dV dV\ dP 1 d{Rr RR ) dr RZ 

p\ VU VV = 1 — — H £i£. (7) 

P \dt dZ OR) OR R OR 8Z y> 

where, r l3 = = //(|^ + ^) , (8) 

f t*o ■ for n<n , 

^ \ an™- 1 + ron- 1 : for U > U K ' 

IT = pE l3 E^ (10) 

The following non-dimensional variables will be introduced in the analysis that follows: 

Z R U V r a aZ +1 P cti , H , b 

^=T> r = — ' u = ~ ' v = ~x> 6 = T' P = — t = — , h— — , = — , 
A a c co A /ic n A A a a 

The equation governing the flow of the fluid can now be written in the form 
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( 9V rr dV .r 9 ^ 

R5 U + U OZ +V OR 



)) 



<9Z 



(13) 



where $ 



n-1 



+r 



,9V 



dU, 



,dU 



2 ^W 2 + ( » )2 + ( - )2} + ( ^ +52 - } 



dZ' 



dR 



dZ' 



(14) 



With the use of the long wavelength approximation (i.e for small 8) and the lubrication approach 
(Shapiro et al 1969, Srivastava and Srivastava 1984, Mishra and Rao 2005, Usha and Rao 
1995, Srivastava and Srivastava 1985), the governing equations stated earlier and the boundary 
conditions describing the flow in the laboratory frame of reference in terms of the dimensionless 
variables when f=0 mav be rewritten as 



dP l ^ilil^ + ro) 

8Z R dR 

OP 

°=-8R 



V> = o, 



du 
dr 



\r or J 



dr 



0, r rz = at r = 0; 



I dip 

u = -— — = at r = h 

r or 



The solution of equation (15) subject to the conditions (17) and (18) is given by 

1 



u(r,z,t) 



(Ph - r ) k+1 - (Pr - 7- ) fc+1 , < r < h 



fc+i 



(k + l)P 

where P = -j* k — -. 

If r be the radius of plug flow region, 



(15) 
(16) 

(17) 
(18) 

(19) 



du 

— = atr 

or 



so that r = To/ P. 
T h at r=h, we find h = T h /P. 



t (say), < t < 1 



If T rz 

Thus ? = 

h T h 

Then we find the plug velocity as 

_ (Ph - Tq)^ 1 
(k + l)P 

In order to determine the stream function, we use the boundary conditions 

tpp — at r = 



(20) 



(21) 
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and ip = ifjp at r = r . 
Integrating (19) and (21), the stream function is found in the form 

pk 



1> 



k+1 



~2 {h r ° ] ~ k + 2 



r < r < h 



and ip p 



P k (h - r ) k+1 r 2 



< r < r 



2(k + l) ' 

The instantaneous rate of volume flow through the tube, Q, is given by 



(22) 
(23) 



Q(z,t) = 2 
P k (h-r ) 



pro rh 

/ rudr + 2 / : 

JO Jr () 



rudr + 2 j rudr 

ro 

' - ' ' 1:2 -J- 3k + 2) + 2r h(k + 2) + 2r, 



l (h 2 (k 2 



Now 



(k + l)(k + 2)(k + 3) 

Q(z,t)(k + l)(k + 2)(k + 3) 



1 

n 



(24) 



dp 

dz 



(h - r ) k+1 (h 2 (k 2 + 3k + 2) + 2r h(k + 2) + 2r 2 ) 
Q(z,t)(k + l)(k + 2)(k + 3) 



^+3(1 _ r )fc+i((jfe2 + sk + 2) + 2r(A; + 1) + 2r 2 ) 
The average pressure rise per wave length is calculated as 



Wjf(i)« 



(25) 



(26) 



One may observe that when n = 1, r = and = 0, the equations (19), (22), (24) and 
(25); the expressions exactly reduce to those obtained by Shapiro et al (1969). Our results also 
match with those of Lardner and Shack (1972), when the eccentricity of the elliptical motion of 
cilia tips is set equal to zero in their analysis for a Newtonian fluid flowing on a uniform channel. 
If we put n = 1 and r = 0, equation (18) tallies with that obtained by Gupta and Seshadri 
(1976) for a Newtonian fluid of constant viscosity. 



4. Quantitative Study 

Theoretical estimates of different physical quantities that are of relevance to the physiological 
problem of blood flow in the micro-circulatory system have been obtained on the basis of the 
present study. For this purpose, the following valid data that are valid in the physiological 
range (Guyton and Hall 2006, Srivastava and Srivastava 1984, Fung 1981, Barbee et al 1994) 
have been used: 

a = 10 to 60/im, = 0.1 to 0.9, f = 0.01 to 0.02, Ap = -300 to 50; r = 0.0 to 0.2; 
Q=0 to 2, n= | to 2. The value of k\ that is considered with the non-uniform geometry 
of the micro-vessel, the value of k\ here has been chosen to match with physiological system. 
Thus for arterioles, which are of converging type, the width of the outlet of one wave length 
is 25% less than that of inlet, while in the case of diverging tubes (e.g. venules), the width 
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r 



Fig. 2: Variation of axial velocity in the radial direction at Z=0.5 



u 




Figs. 3 : Aerial view of the velocity distribution at different instants, n — 2/3, k± — 0, Ap = 
0, r = 0.1, <fi = 0.5 for a non-Newtonian fluid of shear thinning type (a) t=0.0 (b) t=0.25 (c) 
t=0.5 (d) t=0.75. 
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of the outlet of one wave length is 25% more than that of inlet. To obtain the quantitative 
results, the instantaneous rate of volume flow Q(z, t) has been considered to be periodic in (z-t) 
(Srivastava and Srivastava 1984, Gupta and Seshadri 1976, Srivastava and Srivastava 1983), so 
that Q appearing in equation (13) can be represented as 

{Q n + h 2 — 1 — 2 /2 : for sinusoidal wave, 
Q n + h 2 -l- A c 2 /2 - A\ c <p/n : for SSD expansion wavifl) 
Q n + h 2 — 1 — A c 2 /2 + 4A c 0/7r : for SSD contraction wave; 

Q being the time-averaged flow flux. Since the right hand side of equation (13) cannot be 
integrated in closed form, for non- uniform/uniform geometry, for further investigation of the 
problem under consideration, we had to resort to the use of software Mathematica. This helped 
us to calculate the numerical estimate of the pressure rise, Ap. 

4-1. Distribution of Velocity 

For different values of the amplitude ratio 0, flow index number n, k\ and r, Figs. 2-6 present 
the distribution of axial velocity in the cases of free pumping, pumping and co-pumping zones. 
It has been shown in Fig. 2 that the results computed on the basis of our study for the particular 
case of a Newtonian fluid tally well with the results reported by Shapiro et al (1969) when the 
amplitude ratio <fi = 0.5. Since the velocity profiles along with the radius of the tube change 
with time, we have investigated the distribution of velocity at a time an interval of a quarter of 
a wave period. Figs. 3 gives us the aerial view of a few typical axial velocity distribution for 
a Newtonian fluid and a shear-thinning fluid (n=2/3) flowing in uniform micro-vessels, while 
Fig. 4 (cf. Fig. 6 of Takabatake and Ayukawa 1982) shows the velocity distribution in a vessel. 
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Figs. 5 reveals that at any instant of time, there exists a retrograde flow region. However, the 
forward flow region is here predominant since the time averaged flow rate is positive. For a 
shear-thinning fluid (n=2/3), the present study indicates that there exist two stagnation points 
on the axis. For example, at time t=0.25, one of the stagnation points lies between Z=0.0 
and Z=0.25, while the other lies between Z=0.75 and Z=1.0. Similar observations were made 
numerically by Takabatake and Ayukawa (1982) for a Newtonian fluid. Figs. 5(a)-(b) depict 
that in both the regions, as r increases, the magnitude of velocity decreases for both types of 
fluids. It can be observed from Fig. 5(c) that for a Herschel-Bulkley fluid (n=2/3) with Ap = 0, 
the velocity in both the regions of backward and forward flows is enhanced as the value of <f> 
increases when < <fi < 0.6, but beyond <fi = 0.6, the velocity in the backward flow region 
decreases with cf) increasing. 

It is worthwhile to observe the significant influence of the rheological fluid index 'n' on the 
velocity distribution (cf. Figs 5(d)-(i)) for flows in uniform/non-uniform vessels. Figs. 5(d)-(i) 
reveal that the parabolic nature of the velocity profiles is disturbed due to non- Newtonian effect. 
Magnitude of the velocity decreases at the maximum occlusion region; while in the remaining 
regions including backward flow region, a reverse trend is noticed with an increase in the value 
of 'n' (cf Figs. 5(d-f, h-i)). However, in a diverging tube, flow reversal is seen to vanish for 
shear thickening fluids satisfying the rheological fluid index 'n'> 3/2 (cf. Fig. 5(g)). The last 
two observations (cf. Figs. (c,g)) are in contrast to the case of two dimensional channel flow 
was. For a converging vessel, the magnitude of the velocity is greater than that of a uniform 
vessel; however, for a diverging vessel, our observation is altogether different. It is essential 
to mention that Figs. 5(j)-(l) show the influence of pressure on velocity distribution for shear 
thinning /shear thickening fluids. In the case of a uniform/diverging vessel, if Ap decreases, the 
flow reversal tends to decrease; however for a converging tube, although there is reduction in the 
region of flow reversal, the change is not significant for shear thickening fluid. It is important 
to mention that flow reversal due to change in sign of the vorticity or the shear stress along the 
wavy wall. 

The results corresponding to SSD wave is shown in Figs. 5(m-n). It is seen that for SSD 
expansion waves, there is no backward flow and its magnitude is considerable less as compared 
to the sinusoidal wave (cf. Fig. 4). Moreover, as A c increases, velocity is seen to rise except near 
around Z=0.25 where it decreases if A c exceeds 0.6. While for SSD contraction waves, Fig. 5(n) 
shows that if A c = 0.2, there is no backward flow within one wave length. If A c > 0.3, backward 
flow is observed from Z=0.15 to Z=0.7. 

4-2. Pumping Behaviour 

Figs. 6 give the variation of volumetric flow rate of peristaltic waves in the case of our present 
study for different values of the amplitude ratio 0, flow index number n, r, k±. Shapiro et al. [?] 
used lubrication theory to show that the flow rate averaged over one wave varies linearly with 
the pressure difference. For our study of the peristaltic transport of a non-Newtonian fluid, the 
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relationship between the pressure difference and the mean flow rate is found to be non-linear 
(cf. Figs 6). Figs. 6 show that for a non-Newtonian fluid, the mean flow rate, Q increases 
as Ap decreases. Plots in Figs. 6. (a) and (b) indicate that the pumping region enhances with 
an increase in the value of the amplitude ratio, <fi for both shear thinning and shear thickening 
fluids. It has been illustrated in Figs. 6(c-e) that the influence of the rheo logical parameter 'n' 
on the pumping performance in uniform/diverging/converging vessels that the pumping region 
(Ap > 0) significantly enhances with the increase in the value of 'n', while in the co-pumping 
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region (Ap < 0) the pressure-rise decreases when Q exceeds a certain value. As seen in Figs. 
6(f-g), Q is not significantly affected by the value of the parameter r for free pumping case. We 
further find that for both shear-thinning and shear thickening fluids, pumping region increases 
with r increasing. Moreover, when Q exceeds a certain critical limit, Ap = 0, pressure rise 
decreases with the increase in r. 

The effect of SSD wave on pumping is revealed in Figs. 6(h-i). Unlike sinusoidal wave 
form Q increases significantly in pumping, free-pumping as well as co-pumping regions for SSD 
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TO 



Figs. 5: Distribution of velocity in different situations 



expansion wave fronts as A c increases. For SSD contraction wave fronts, our observation is 
altogether different against SSD expansion wave fronts. However, Q < when Ap < for 
shear-thinning fluids in the case of contraction waves. This observation is in contrast to the case 
of sinusoidal was. 

4-3. Streamlines and Trapping 

It is known that one of the important characteristics of peristaltic transport is the phenomenon of 
trapping . It occurs under some conditions as streamlines on the central line are split to enclose a 
bolus of fluid particles circulating along closed streamlines in the wave frame of reference. Then 
the trapped bolus moves with a speed equal to that of the wave. This physical phenomenon may 
have responsibility to the formation of thrombus in blood. Under the purview of the present 
study, Fig. 7-10 give an insight into the changes in the pattern of streamlines and trapping that 
occur due to the changes in the values of different parameters that govern the flow of blood in 
the wave frame of reference. Figs. 7 provide the streamline patterns and trapping of a shear- 
thinning fluid for different values of 0. With an increase in 0, the bolus is found to appear in a 
distinct manner. Streamlines for different fluid indices 'n'are depicted in Figs. 8. These figures 
indicate that occurrence of trapping is strongly influenced by the value of the fluid index. Figs. 
9 show that trapped bolus increases in size and also that it has a tendency to move towards the 
boundary as the flow rate increases. It is important to note that the bolus appearing for r = 
is going to disappears as the value of r rises (Figs. 10). 



4-4- Distribution of Wall Shear Stress 

It is known that if the shear stress generated on the wall of a blood vessel exceeds a certain 
limit, the constituents of blood are likely to be damaged. Also, the magnitude of the wall 
shear stress has a vital role in the process of molecular convection at high Prandtl or Schmidt 
number (Higdon 1985). In view of the above observations, it is important to study the shear 
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Figs. 6: Pressure rise versus flow rate 
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Figs. 7 : Streamline patterns for peristaltic flow of a shear-thinning fluid for different values of 
(n=2/3, Q = 0.4, h = 0, r = 0) 
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Figs. 8 : Streamline patterns and trapping in the case of peristaltic flow for different values of 
physiological fluid index V when = 0.5, r = 0.1, k-y — 0, Q — 0.4 



stress that is developed during the hemodynamical flow of blood in arteries. The wall shear 
stress distributions are plotted in Figs 11 under varied conditions. The distributions of wall 
shear stress at four time instants during one complete wave period have been presented in Fig. 
11(a). It may be observed from this figure that at each of these instants of time, there exist two 
peaks in the wall shear stress distribution, with a gradual ramp in between; however, negative 
peak of wall shear stress, T min is not as large as the maximum wall shear stress, r max . The 
transition from T min to r max of the wall shear takes place in some zone between the minimum 
and maximum radii of the vessel. At the location where the maximum occlusion occurs, wall 
shear stress along with the pressure is maximum. Since the pressure gradient to the left of this 
location takes positive value, the local instantaneous flow will take place towards the left of 
Tmax- This may lead to some consequences may be happen, for example, if the shear rate at the 
crest exceeds some limit, a dissolving wavy wall will have a tendency to level out. Moreover, 
some bio-chemical reaction between the wall material and the constituents of blood may set in. 
As a result of this, there may be deposition of the products of the chemical reaction on the 
endothelium and consequently wall amplitude may increase at a rapid rate. This is likely to 
lead to clogging of the blood vessel. 
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Figs. 9 : Streamline patterns and trapping for the effect of Q when n — 2/3, 
0,r = 0.1 
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The peaks of the wall shear stress distribution on both sides of r max are small and hence 
the local instantaneous flow will occur in the direction of the peristaltic wave. For a Herschel- 
Bulkley fluid, Figs, ll(b-c) show that in the contracting region where occlusion takes place, 
there is a remarkable increase in the wall shear stress due to an increase in the value of <fi for 
shear-thinning as well as shear-thickening fluids. With the increase in 0, magnitude of r min 
increases in the expanding region for shear-thinning/ shear-thickening fluids, although the effect 
is not very prominent. It may be observed from Fig. 11(d) that r max increases with increase in 
r, while r has little effect upon r min . 

The quantum of influence of the rheological fluid index 'n' on the distribution of wall shear 
stress is shown in Figs, ll(e-g) for uniform/non- uniform blood vessels. In all types of vessels 
studied here, r max increases with the increase in 'n'. Moreover, it is very important to mention 
that the shear stress difference between the outlet and the inlet in the case of a converging vessel 
is exceedingly large in comparison to the case of a diverging vessel. As the time averaged flow 
rate increases, Figs, ll(h-j) indicate very clearly that the wall shear stress tends to decrease for 
all the types of vessel examined here. One can observe from Figs, ll(k-l) that r h changes its 
values within the region of SSD wave activation; beyond this , however it remains constant. 
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Figs. 10: Streamline patterns for different values of r in the case of peristaltic flow when 
n = 2/3, = 0.5, h =0, g = o.4 



5. Summary and Conclusion 

The present paper deals with the peristaltic motion of blood in the micro-circulatory system, by 
taking into account the non-Newtonian nature of blood and the non-uniform geometry of the 
micro-vessels, e.g. arterioles and venules. The non-Newtonian viscosity of blood is considered 
to be of Herschel-Bulkley type. The effects of amplitude ratio, mean pressure gradient, yield 
stress and the rheological fluid index n on distribution of the velocity and wall shear stress, 
pumping phenomena, streamline pattern and trapping are investigated under the purview of the 
lubrication theory. Experimental observations have revealed that in the case of roller pumps, the 
fluid elements are quite prone to significant damage; also during the process of transportation 
of fluids in living structures executed by using arthro-pumps, the fluid particles are likely to 
be appreciably damaged. Qualitative and quantitative studies of the current analysis for the 
wall shear stresses have a significant bearing in extra-corporal circulation, where the heart- 
lung machine is usually used and thereby there is a chance of the erythrocytes of blood being 
damaged. Thus the present study bears the potential of significant application in biomedical 
engineering and technology. 
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The present study reveals that at any instant of time, there is a retrograde flow region 
for Herschel-Bulkley type of non-Newtonian fluids like blood when Ap < 0. The regions of 
forward/retrograde flow advance at a faster rate, if the values of n and <fi are raised. In the case 
of a uniform/diverging tube, the flow reversal tends to decrease when Ap tends to be negative; 
however, for a converging tube such an observation is not very prominent. Moreover, the present 
study shows that non-uniform geometry affects quite significantly the distributions of velocity 
and the wall shear stress as well as pumping and other flow characteristics. The amplitude ratio <fi 
and the rheological index 'n' are very sensitive parameters in changing the peristaltic pumping 
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Figs. 11: Wall shear stress distribution in different situations. 



characteristics and the distribution of velocity and wall shear tress. Moreover, the parabolic 
nature of the velocity profiles is found significantly disturbed by rheological fluid index V. 

From the present study we may conclude that for SSD expansion waves, the pumping per- 
formance is better than that in the case of sinusoidal waves. It is also important to note that 
the backward flow region is totally absent in the case of SSD expansion waves. Further, it may 
be concluded that backward flow originates from the contraction of vessels. 
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